Bulk Mediated Surface Diffusion: Non Markovian Desorption 

with Finite First Moment 

Jorge A. Revelli, 1 ]*! Carlos. E. Budde, 1 !] Domingo Prato/'P and Horacio S. Wio 2 ' 3 ^ 

1 Facultad de Matemdticas, Astronomia y Fisica, Universidad Nacional de Cordoba 

5000 Cordoba, Argentina; 
2 Instituto de Fisica de Cantabria, E- 3 9005 Santander, Spain 
3 Grupo de Fisica Estadistica, Centro Atomico Bariloche and Instituto Balseiro, 
84 00 San Carlos de Bariloche, Argentina 

Abstract 

Here we address a fundamental issue in surface physics: the dynamics of adsorbed molecules. 
We study this problem when the particle's desorption is characterized by a non Markovian process, 
while the particle's adsorption and its motion in the bulk are governed by a Markovian dynamics. 
We study the diffusion of particles in a semi-infinite cubic lattice, and focus on the effective diffusion 
process at the interface z = 1. We calculate analytically the conditional probability to find the 
particle on the z = 1 plane as well as the surface dispersion as functions of time. The comparison 
of these results with Monte Carlo simulations show an excellent agreement. 
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I. INTRODUCTION 



The mechanism called bulk-mediated surface diffusion typically arises at interfaces sepa- 
rating a liquid bulk phase and a second phase which may be either solid, liquid, or gaseous. 
Whenever the adsorbed species is soluble in the liquid bulk, adsorption-desorption pro- 
cesses occur continuously These processes generate surface displacement because desorbed 
molecules undergo Fickian diffusion in the liquid's bulk, and are latter re-adsorbed elsewhere. 
When this process is repeated many times, an effective diffusion results for the molecules 
on the surface. The importance of bulk-surface exchange in relaxing homogeneous sur- 
face density perturbations is experimentally well established . Adsorption at 
the solid-liquid interfaces arises, for instance, in many biological contexts involving protein 




|8(, in solutions or melts of synthetic macromolecules 



ff 



121, in 



deposition 

colloidal dispersions Q], and in the manufacture of self-assembly mono- and multi-layers 

Usually the studies performed in this type of systems are done within the framework of 
a Master Equation scheme jisl. tlflj|. where the particle's motion through the bulk and the 
adsorption-desorption processes are Markovian. In two recent papers |2CJ,|2jJ we have shown 
the most important features of these phenomena. Non Markovian diffusion processes have 
been used in order to modelling a great diversity of phenomena some of the most relevant 
are: anomalous charge transport in amorphous materials (3], diffusion of particles with 
internal states j^, chromatography [3], dielectric relaxation due to defect diffusion j^j]. 

We address one of the fundamental issues in surface physics: the dynamics of adsorbed 
molecules. We study this problem when the particle's desorption is characterized by a non 
Markovian process, while the particle's adsorption and its motion in the bulk are governed 
by a Markovian dynamics. We analyze the diffusion of particles in a semi-infinite cubic 
lattice, and focus on the effective diffusion process at the interface z — 1. We calculate 
analytically the conditional probability to find the particle on the z = 1 plane as well as the 
surface dispersion as functions of time, and test those results comparing with Monte Carlo 
simulations. When non Markovian processes are present Generalize Master Equations are 
involved, these equations being characterized by a "memory kernel" that is related univocally 
with a Continuous Time Random Walk (CTRW) scheme j^ . 

Non Markovian desorption has been introduced in many fields of natural science de- 
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scribing transport processes in chemical reactions [2q . re-emission when surfaces contains 
"deep traps" |27(, capture and re-emission from surfaces that contain sites with many inter- 
nal states such the "ladder trapping model" proteins with active sites deep inside its 
matrix , etc. 

The main goal of this work is to get information about the influence of a non Markovian 
desorption dynamics on the effective diffusion process at the interface z = 1, and in this way 
to develop some criteria for looking for non Markovian desorption effects in experimental 
situations. For that purpose we calculate analytically the temporal evolution of the variance 
((r 2 (t)) plane) and the conditional probability of being on the surface at time t since the 
particle arrived there at t — (P(x, y, z — 1; t\0, 0, 1, t — 0)), that we indicate by P{z = l,t). 

When ip(t), the waiting time density for desorption defined in the CTRW scheme, has a 
finite first moment we have been able to establish the long time asymptotic behavior for 
(r 2 '(t)) pi ane and P(z = l,t) for any non-Markovian desorption. This behavior is the same 
as the Markovian one and only depends on the first moment of ip{t). We also analyzed 
the behavior associated to two different families of desorption waiting time densities. For 
the first one we have observed two regions, called strong and weak adsorption regime (to 
be defined later). In the strong adsorption limit we observe a transient regime, where the 
temporal evolution is characterized by damped oscillations (its frequency not depending on 
the degree of departure from Markovian behavior), and a Markovian asymptotic one. The 
oscillatory effects disappear in the weak adsorption region. For the second kind of desorption 
dynamics, a transient regime with a non-monotonic behavior for the slope of the variance 
emerges, the asymptotic behavior being again Markovian like. 

The organization of the paper is as follows. In the next Section we formally present the 
model, supported by a Generalized Master Equation which describe the particle's dynamics 
through the bulk and the adsorption-desorption dynamics at the surface. We devote Section 
III to present and discuss the time dependence of (r 2 (t)) p ; ane and P{z = 1, t) for two families 
of specific memory kernels, and compare the analytical results with Monte Carlo simulations 
in the first case, and discuss the general characteristics in the second. In Section IV we 
study the asymptotic behavior of the aforementioned magnitudes, in the asymptotic long 
time limit. In the last Section we discuss the results and present some conclusions. 
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II. THE ADSORPTION-DESORPTION MODEL 



Let us start with the problem of a particle making a random walk in a semi-infinite cubic 
lattice (with a lattice constant equal to one). The position of the walker is defined by a 
vector r whose components are denoted by a set of integer numbers n, m, I corresponding 
to the directions x, y and z respectively. The probability that the walker is at (n, m, I) for 
time t given it was at (0, 0, l ) at t — 0, P(n, m, I; t\0, 0, l , t = 0) = P(n, m, I; t), satisfies the 
following Generalized Master Equation 

P(n,m,l;t) = jP(n,m,2;t) - [ dt' K(t')P(n,m,l;t - t'), for I = 1 

Jo 

P(n,m,2;t) = / dt' K(t')P(n,m,l;t - t') + ^P(n,m,3;t) - ^P(n,m,2;t) 
Jo 

+a[P(n - 1, m, 2; t) + P(n + 1, m, 2; t) - 2P(n, m, 2; t)} 
+P[P(n, m - 1, 2; t) + P(n, m + 1, 2; t) - 2P(n, m, 2; t)], for I = 2 
P(n, m, I; t) = a[P(n -l,m,l;t) + P(n + 1, m, I; t) - 2P(n, m, I; t)) 
+P[P(n, m-l,l;t) + P(n, m + 1, 1; t) - 2P(n, m, I; t)) 
+^[P(n,m,l + l;t) + P(n,m,l-l;t)-2P{n,m,l;t)}, for I > 3, (1) 

a, (3 and 7 are the transition probabilities per unit time through the bulk in the x, y and z 
directions respectively and K(t) represents the memory kernel at all sites (n, m, I = 1) and 
(n, m, I — 2) . The form of these equations is similar to those indicated in |2£], |21|, with the 
desorption parameter S replaced by the kernel K(t). 

In order to solve the above equations we follow the same procedure as in [3, 12]] , taking 
the Fourier transform with respect to the x and y variables and the Laplace transform with 
respect to the time t in the above equations, we obtain 

sG(k x , k y , 1; s) - P(k x , k y , 1, t = 0) = ^G(k x , k y , 2; s) - K(s)G(k x , k y , 1; s), 1 = 1 
sG{k X i ky, 2\ s) P[k X: ky, 2, t = 0) = A{k x , ky)G{k x , k y , 2;s)+K(s)G(k x ,k y ,l;s) + 

-fG(k x , k y , 3; s) - 2^G{k x , k y , 2;s), 1 = 2 
sG(k x , ky, l] s) P(k x , ky, l,t = 0) = A(k x , k y )G( y k x , ky, /; s) + 

j[G(k x , k y , I - l;s) + 

G(fc a>J Ai / ,Z + l;s)-2G(A; iB ,Ai /J Z;a)]. />3 (2) 
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We have used the following definitions [20|, |21 1 



G(k x ,k y ,l;s) = G(k x ,k y ,l;s\0,0,l ;t = 0) 



e -st e i( - k * n+k y m) P(n,m,l;t)dt 

n,m,— 00 

00 

£[ ^ e^ n+feam) P(n,m,/;t)], 



(3) 



n,m,—oo 



where L indicates the Laplace transform of the quantity within the brackets, and |20L \21\ 



A(k x , k y ) = 2a[cos(k x ) - 1] + 2/3[cos(/c 3/ ) - 1]. 
Equation (0) may be expressed in matrix form as 

[si - H}G = I, 
where the square matrix G has components 

G ll0 = [G[k x ,ky,l;s\0,0,l o ;t = 0]]. 
I is the identity matrix and if is a three-diagonal matrix with the following form 

j -K(s) 7 



(4) 



(5) 



(6) 



H 



K{s) C 7 •■• 
7 C 7 •■• 
7 C 7 ••• 

y 



where C is defined as 



C = -2 1 + A(k x ,k y ). (7) 

n 

Above equations for G(k x , k y ,l; s) are similar to Eq.(2) and Eq.(5) in [2(j] with K(s) 
replacing 5. Therefore, all results obtained in that paper remains valid for non-Markovian 
dynamic when 5 is replaced by K(s) in the Laplace domain. 

The expression for the Laplace transform of the variance ((r 2 (s)) p i ane ) became 

N(s) 



L[(r (t)} p l ane ] = (f (s)) p l ane 



D(s) 



(8) 
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with 



N(s)= [4K(s)j(a + /3)} 

2 7 3 - s 2 ( v / s(47 + s) - s) - 37 2 (v / s(4 7 + s) - 3m) - 2 7 s(2 v/s(4 7 + s) - 3s) 



(9) 



D(s) = v/s(4 7 + s) 



(10) 



7s(2 7 + s - y/s(A^ + s)) + fsT(s) (7(^5(47 + s) - 3s) + s(^/s(4-? + s) - s)) 



The Laplace transform of P(z = 1, s), the conditional probability of being in the surface, 



is 



L[P(z = 1, t)] = P(z = l,s)= 7(27 + s - VK47 + s 



(11) 



7s(2 7 + s - y/s(4ry + s)) + K(s)(j(yG{^Ts) - 3s) 
+s(v/s(4 7 + s) -s)) . 

Finally we want to point out that the relation between ip(t), the waiting time density for 
he desorption processes defined in the CTRW scheme, and the memory kernel of Eq. ((TJ 
25I ] in the Laplace domain is given by 

K(s) 



;i-#o)" 



(12) 



III. ANALYTICAL RESULTS AND MONTE CARLO SIMULATIONS 



Here we show the results obtained analytically and compare them with Monte Carlo 
simulations. In all cases we have fixed the parameters: a — [3 — 7 = 1. The simulations 
results were obtained averaging over 2 x 10 6 realizations. 

To describe the desorption dynamics from the surface we have chosen two families of 
waiting time densities (ip{t)). The first of them was introduced by Scher and Lax 2^] to 
describe the frequency dependence of the electric conductivity in disordered solids when 
transport is due to impurity hopping. It has been extensively exploited in modelling non 
Markovian cases that emerge when an average over transition rates (disorder) is taking 



into account 



22] and 



The reason of its wide use are its versatile functional form and 



its simplicity which allows to take into account a controllable spread of transition rates 



22|. When only one transition rate is present the Markovian description is reobtained (the 



memory kernel is a Dirac delta function) and when the spread is very wide, steps occur at 
fixed regular intervals of time (see below). The adopted function (see also 30j) is 

m = 9a^^e^\ (13) 

where a is a positive integer and T(a) is the Gamma or Factorial Function. It is worth 
remarking here two important facts about this family of functions. First, as can be seen 
from the Eq. 1)1 there are two parameters which characterized the function. The param- 
eter a, called markovianicity parameter, defines the degree of function's departure from the 
Markovian behavior (a = 1 corresponds to the Markovian case; a ^ 1 to the non-Markovian 

n 

case), while the parameter 9 is the "average desorption's rate". Second, as shown in [30], 
the mean value of these waiting time densities is 

poo 

(t) = t V(0 dt = 0"\ (14) 
Jo 

showing that the "average desorption's time" does not depends on the parameter a. For the 
form of this family of functions, see Fig. 1 in [30]. 



In Figs. QandElwe present the temporal evolution for the conditional probability (P(z = 
l,t)) for the Markovian and two non Markovian cases (two different values of a) with the 
same value of 9. From these figures we can observe two temporal regions: a transient one 
which ranges from t = to t « 1000 and an asymptotic region (t > 1000) where the behavior 
approaches that of the Markovian case. It is important to remark the existence of damped 
oscillations in the transient region for the non Markovian case. This oscillatory behavior is 
due to the waiting time density functions used as will be explained later. 

In Figs. 01andE]we depict the results for the variance (r 2 (t)) p i ane , for both the Markovian 
and the two different non Markovian cases (a = 75 and a = 100). We can observe two 
important features. The first one is that the variance in the non Markovian regime shows 
a delay in the beginning of the spreading. The second feature is that the system presents 
an oscillatory-like behavior, with the oscillations attenuating as the time grows. This fact 
is still more apparent in the insert of the figures where we show the temporal evolution of 
the spreading velocity {Vspread)- It is worth remarking here that oscillations only appears in 
the non Markovian case with a ^> 1 and are due to the particular behavior of the family 
of waiting time density defined by Eq. (|13j). As it is known, this function goes to a Dirac 
delta function as the markovianicity parameter a —>■ oo (ip(t) =^ 5(t — 0~ 1 )), implying a 




FIG. 1: Temporal evolution of the P(z = l,t). We have represented the case for 8 = 0.01. Dot 
line depicts Markovian evolution, meanwhile continuous line and open circles (which correspond 
to Monte Carlo simulations) represent the non Markovian case for a = 20 

periodic-like behavior. This "periodicity" explains the oscillatory behavior of P(z = l,t) 
and {r 2 (t)) p i ane . Also the figures show an excellent agreement between the theoretical and 
Monte Carlo simulation results. 

We remark here that we can distinguish two region defined by the ratio 9/j. The strong 
adsorption region characterized by ^ << 1 where the non Markovian character shows im- 
portant differences respect to the Markovian case in the transient temporal regime and the 
weak adsorption region > 1) where these differences disappear. 

Another worth remarking point from these figures is that u, the frequency in the oscil- 
lation, remains unperturbed due to the fact that all desorption waiting time densities used 
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II 
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FIG. 2: Temporal evolution of the P(z = l,t). We have represented the case for 9 = 0.01. Dot 
line depicts Markovian evolution, meanwhile continuous line and open circles (which correspond 
to Monte Carlo simulations) represent the non Markovian case for a = 75 



have the same "average desorption's time". This aspect becomes apparent in Fig. |S| 

The second class of desorption dynamics that we consider in this work corresponds to 
the so called "direct-access multi-trapping process" [28]. Such an approach was proposed 
in order to describe in a phcnomenological manner the transport of excited charge carriers 

nn n 

across amorphous material |3JJ,|32J], or the motion of intersticial in metals [33], etc. When 
the walker arrives at a surface site this model considers transitions to and from (N — 1) 
internal states, both with probability A per unit time, before it desorbs to the bulk with 
probability 5. The resulting desorption probability waiting time density in Laplace domain 
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FIG. 3: Temporal evolution of the (r 2 (t)} p i ane . We have taken 9 = 0.01. The solid line and open 
circles (which correspond to Monte Carlo simulations) depict the non Markovian case with a = 75. 
In the insert we can see the Vs pre ad vs. t for this case (solid line) and the Markovian evolution (dot 
line) 



is 

^ = i^^m (l5) 

and the "average desorption's rate" in this case is (t) = y. For a detailed analysis of this 
model and derivation of Eq. ()15|) see Ref. {^J. In Figs. El and [7| we depict the results for the 
variance (r 2 (t)) p i ane for two different non Markovian cases (N = 2 and N = 10 respectively). 
We also show the Markovian behavior with the same average time for each case ((£) =2 
for N = 2 and (t) = 10 for iV = 10). We can observe that the temporal evolution of the 
variance in the non Markovian case shows a transient regime characterized by the slope's 
non-monotonic behavior, which increases with N. This fact is still more apparent in the 
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FIG. 4: Temporal evolution of the variance (r 2 (t)) p i ane . We have represented two cases. Dot line 
represents the Markovian case and the continuous line depict the non Markovian (a = 100) for 
6 = 0.01. In the insert the curve denoted by a represents the Markovian case; the b curve represents 
the non Markovian one. 

insert of the figures where we show the temporal evolution of the spreading velocity (Vspread)- 

IV. ASYMPTOTIC BEHAVIOR 

The asymptotic behavior for large t of (r 2 (t)) p i ane and P{z = l,t) can be obtained 
analyzing the limit s C 1. We assume that ip(t), the waiting time density for desorption, has 
a finite first moment (t) which means that ip(t) ~ 1 — (t)s when s < 1 and consequently 
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FIG. 5: u vs. 0. Open circles represents a = 100, squares a = 200 and cruxes a = 500. 



K(s) ~ in this limit. From Eqs. f 

(r 2 (t)) pla 



ij, (HOI) and (PJ we obtain 

V7(« + /?)(t) 



P(* = l,t) 



2s« 



(16) 



(17) 



,S'2 



Exploiting Tauberian theorems [25|, the behavior for large t can be obtained resulting in 



(r 2 (t)) pl0 



2r[3/2] 



*3. 



v ' ; r[i/2] 



(18) 



(19) 



Equations (fTHj) and (fH?j) represent the behavior of (r 2 (t)) p i ane and P(z = for any non- 
Markovian desorption for large values of time. This behavior is the same as the Markovian 
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FIG. 6: Temporal evolution of the (r 2 (t)) p i ane . The solid line depicts the non Markovian case with 
N = 2, 5 = 1 and A = 0.01 and the dot line the Markovian one. In the insert we can see the 
y Spread vs - t for this case (solid line) and the Markovian evolution (dot line) 

one and only depends on the first moment of ip(t). It is important to remark that the only 
assumption was that ip(t), the waiting time density for desorption, has a finite first moment 
(t). 

V. CONCLUSIONS 

We have studied here the evolution of particles diffusing on a surface. The diffusion have 
been performed across the bulk surrounding the surface, this phenomenon being called bulk 



mediated surface diffusion\2(t l21| . Usually the proposed models are based on Markovian 
desorption processes. The main feature of this work was to present an analytical model 
for non Markovian desorption from the surface. The bulk that surrounds this surface was 
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FIG. 7: Temporal evolution of the (r 2 (t)} p i ane . The same parameters as in Fig. 6 but with N = 10 

considered to be semi-infinite and that the particles undergo a Markovian motion on them. 

We observed the influence of the non Markovian desorption dynamic on the effective 
diffusion process at the interface z = 1 by calculating analytically, in the Laplace domain, 
the temporal evolutions of the variance ({r 2 (t)) p i ane ) and the conditional probability of being 
on the surface at time t since the particle arrived there at t — (P(z = l,i)). When the 
waiting time density for desorption has a finite first moment we have established the 
behavior of the above magnitudes for any non-Markovian desorption for large values of time. 
This behavior is the same as the Markovian one and only depends on the first moment of 



We have chosen two families of non-Markovian desorption waiting time densities and 
tested the analytical results for (r 2 (t)) plane and P(z — 1, t) by comparison with Monte Carlo 
simulations, obtaining an excellent agreement. For the first desorption waiting time density 
we can establish two regions based on the ratio calling them strong adsorption (this 
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occurs when 0/7 << 1) and weak adsorption (when 0/7 > 1). In the strong adsorption 
limit we can observe a transient regime and a stationary one. The main feature in the 
transient regime is that the temporal evolution is characterized by damped oscillations whose 
frequencies are in direct relation to 6, the average desorption rate. It is worth remarking 
here that the frequency does not depend on the markovianicity parameter. The effect of 
this parameter appears on the amplitude of the oscillations which disappear in the weak 
adsorption region. For the second kind of desorption dynamics we found a transient regime 
with an emerging non-monotonic behavior for the slope's variance. In the asymptotic regime, 
the conditional probability of the non Markovian system tends to the Markovian one as 
is expected from the analytical results obtained in Section IV. In a recent paper |24j we 
have shown that, when we consider finite or infinite biased systems, and desorption 
waiting time densities with finite first moment, the variance growths linearly with time in 
the asymptotic regime. 

Finally, it is worth here remarking an important aspect of the present approach. We have 
shown through the above results that the behavior of (r 2 (t)) p i ane and P(z = l,t) results to 
be strongly dependent on the desorption mechanism. As the effective dispersion and the 
percentage of particles that remain on the plane z = 1 are measurable magnitudes they 
could be used to investigate the characteristic and get information about the fundamental 
parameters of the desorption processes. 
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